Gaussian Statistics of Fracture Surfaces 
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We analyse the statistical distribution function for the height fluctuations of brittle fracture 
surfaces using extensive experimental data sampled on widely different materials and geometries. 
We compare a direct measurement of the distribution to a new analysis based on the structure 
functions. For length scales 5 larger than a characteristic scale 8* , we find that the distribution of 
the height increments Ah — h(x + S) — h(x) is Gaussian. Self-affinity enters through the scaling of 
the standard deviation a, which is proportional to 5^ with a unique roughness exponent. Below the 
scale 5* we observe an effective multi-affine behavior of the height fluctuations and a deviation from 
a Gaussian distribution which is related to the discreteness of the measurement or of the material. 

PACS numbers: 83.80.Ab, 62.20.Mk, 81.40.Np 



It is difficult to believe that there may be anything 
in common between the morphology of fractures in, say, 
concrete and aluminium, except for the qualitative state- 
ment that they both are "rough". The roughness seems 
very different when comparing the two materials. Stud- 
ies 0, 0, 0, 0, Hi have shown that the scaling properties 
of this roughness are the same to within the measuring 
accuracy for not only these two materials, but for most 
brittle or weakly ductile materials that have been tested. 
The scaling properties of the roughness alluded to above, 
is more precisely described as the fractures being self- 
affine. The typical deviations Ah of the surface as a 
function of distance S along the fracture surface scale as 
Ah<x5 c [10H 0,0. K has been suggested that these 
scaling properties might be universal 

Most studies focus only on the scaling properties of the 
fracture surfaces. They give no hint of the actual statisti- 
cal distribution giving rise to such a scaling. In this study 
we go beyond a calculation of the roughness exponents 
and propose a statistical distribution for the height fluc- 
tuations Ah of fracture surfaces. For the various materi- 
als and geometries we analyzed, we find that the Gaus- 
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sian distribution provides a complete statistical descrip- 
tion of the morphology of fractures at least at large scales. 
This result is in contradiction with other works where 
multi-scaling_J| , more precisely, global multi-affinity is 
observed 0,H. 

Our work is based on the analysis of experimental data 
obtained from various experiments on different materials. 
The materials have been broken in different modes and 
geometries and the surfaces have been analyzed along one 
or more directions. 

First we used the roughness measurement of a fracture 
surface obtained from the failure of a granite block in 
mode I (4 bending point failure). The scanned area of 
lOcmx 10cm covers the complete section of the block with 
a grid mesh of (<5 ) 2 =48/mix48/im. Accordingly the grid 
size is: 2062x2063, i.e. more than 4 million data points. 
The profiler is optical with a laser beam of 30/im in diam- 
eter |9(. To reduce possible optical artifacts, the scanned 
surface comes from a high-resolution silicon mold of the 
granite fracture. The replica technique in a perfectly ho- 
mogeneous material removes fluctuations of local optical 
properties and significantly improves the quality of the 
roughness measurement. 

Second, we used data from interfacial fracture fronts 
propagating (in mode I) into the annealing plane of two 
plexiglas (PMMA) plates [H IH H 13 • We have ana- 
lyzed 6 long front lines (obtained by assembling fronts for 
a crack at rest ^(|) containing 17000 pixels each, with 
a pixel size 8 = 2.6/im. The roughness exponent was 
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FIG. 1: Convergence of the moment ratios Rt as function FIG. 2: Data collapse of the structure function normalized 



of S towards the Gaussian ratios R^ for the different order 
k = 1, 2, 6 , confirming the underlying Gaussian statistical 
distribution. For each set, on the small scales, the individ- 
ual lines represents from below to above increasing k values. 
The various ratios Rk are averaged over a) 6 interfacial fronts 
in PMMA, b) 5 fracture fronts in Paper and 100 profiles c) 
perpendicular to the fracture propagation in a granite block. 



by the Gaussian ratios Ck(S)/R^ for the values k = 1, 2, 6. 
The various data set are displaced vertically to improve the 
visual clarity. The dashed lines are fit to the second order 
structure functions C2(S)/ R§ on a range where the structure 
functions collapse. Their slopes provide an estimate of the 
roughness exponent. We estimate the following exponents, for 

■ 0.8, 



PMMA C « 0.5, for paper C 2d 
respectively. 



0.6 and for granite £ 



found equal to C = 0.63 ± 0.03 (TfJ- 

Finally, we studied fracture surfaces obtained during 
fracture experiments on a quasi two-dimensional mate- 
rial: fax pap er sheets loaded in mode I at a constant 
force fLU Il5j . High resolution scans were performed on 
post-mortem samples. We analyzed 5 fronts with around 
10000 pixels each, the pixel size is S a = 20/im. In that 
case the roughness exponent measured is ( 2d « 0.6. 

We aim at estimating the statistical distribution func- 
tion of the height fluctuations of a fracture surface 
P(h(x + S) — h(x)). However, a direct measurement of 
the distribution function is not always accessible due to 
limited statistics. Only the first data set which is very 
large, will allow us to perform such a direct estimate. For 
the others, we propose a method introduced in con- 
nection with the study of directed polymers that is 
based on the structure functions defined as the k th root of 
the k th moment of the increment \Ah\ = \h(x + S) — h(x)\ 
on a scale 6 : 



C k {5) = {\h(x + 5)-h{x)\ k ) 1 ' k 



(1) 



The average is taken over the spatial coordinate x. Now, 
forming the ratio between the k th structure function and 
the second structure function, we define the function 



Rk(S) = 



(\h(x + s) - h{x)\ k y/ k 

((h(x + 5) - h(x)) 2 } 1 / 2 



(2) 



In the case of a one-dimensional Brownian motion, the 
statistical distribution of increments Ah is a Gaussian: 



P(Ah) = 



1 
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(3) 



Note that the self-affinity enters through the variance 
of the distribution a 2 oc S 21 *, where £ = 1/2 for the 



Brownian motion. Below we find the same statisti- 
cal distribution for the various fracture fronts investi- 
gated except that generally £ ^ 1/2. For the Gaussian 
distribution the moments Eq. Q are easily calculated, 

C°{8) = (25 2 ?) 1 / 2 (T((fc + l)/2)/V^) 1/fe • In this case 



the ratios of the structure functions become 




(4) 



Note that these Gaussian ratios are independent of a 2 
and 5; they contain no adjustable parameters. The gen- 
erality of this result transcends the derivation based on 
an underlying Gaussian distribution. A different under- 
lying distribution will give rise to a different set of ratios 
Rk- 

First, we computed the moment ratios normalized by 
the Gaussian values, Rk{S) / R^ , and in particular their 
variation with S for the different fracture profiles we have 
studied. In Fig.^ we show that the ratios Rk(S) converge 
for different values of k and for large S, towards the values 
of a Gaussian process, suggesting the Gaussian nature of 
the distribution for a range of 5 values. Note that the 
deviations from a constant value for the largest values of 
6 happen when the statistics become very poor. 

Second, we examined the scaling behavior of the struc- 
ture functions. In Fig.[3we show directly that the struc- 
ture functions collapse, when normalized by the Gaussian 
ratios, Ck(S)/R^. The collapse at the larger scales pro- 
vides clear evidence that no multi-scaling is present, i.e. 
the scaling exponent of Ck(S) ~ S^ k is independent of k 
and therefore we may extract a unique roughness expo- 
nent Cfe = Q. A fit to the second order structure functions 
C2(<5)/i?2 on the range where the data collapse provides 
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an estimate of £ within 10%. For the granite surface we 
find a roughness exponent around ( 3d « 0.8. For the in- 
tcrfacial fracture fronts in PMMA, we obtain a roughness 
exponent £ « 0.5, which is slightly lower than but con- 
sistent with previous estimates 0, 0, 0, 0| . Finally, 
for paper we observe an exponent C, 2d ~ 0.6. 

The fact that the rescaling by Gaussian ratios leads to 
a data collapse, suggests that the underlying distribution 
is Gaussian. This result is confirmed by a direct analysis 
of the large data set from the fracture surface in a gran- 
ite block. The analyzed data set consisted of 2000 x 2000 
points representing the central part of a 3D map of the 
fracture surface to reduce boundary effects. From the 
map we not only computed the structure functions as 
shown in Fig. [21 but we also computed directly the sta- 
tistical distribution of the height fluctuations P(A'h) at 
different length scales S, where A'h = (Ah — (Ah))(5). 
Note that we subtract the averaged height fluctuations 
(Ah) in order to center the various distributions around a 
zero mean. The structure functions (Eq. arc defined 
without such a procedure. However, we checked that it 
did not influence the scaling behavior of the structure 
functions, by directly detrending the various profiles. In 
Fig. |3 we show the distributions of the height fluctua- 
tions for logarithmically increasing length scales S. The 
data were extracted in the direction perpendicular to the 
fracture propagation and the distributions were sampled 
from the 2000 profiles h(x) each containing 2000 points. 
We clearly see that above a characteristic length scale 
8* ~ 50 x 48/im the shape of the distributions become 
Gaussian. This scale corresponds to the point of conver- 
gence observed in Fig. ^ and the onset of collapse in Fig. 
121 Interestingly, a similar cross over has been observed in 
the width distribution of contact lines measured recently 
[l8j . We emphasize that the sclf-affine behavior of the 
fracture front enters through the scaling of the standard 
deviation a oc 5< with C 3d « 0.75, sec the inset in Fig. 

m 

Finally, the length-scale independence of Eq. Q and 
the data collapse in Fig. [21 provide clear evidence that 
there is no multi-scaling of the structure functions. That 
being said, we do observe a separation of the structure 
functions at small scales together with a broadening of 
the tails of the distributions observed in F ig, p Sev- 
eral authors have reported similar findings [HJ, [2(| . In 
KPZ models on kinetic surface roughening a broadly dis- 
tributed noise gives rise to rare but large perturbations 
of the surface and hence a multi-scaling of the structure 
functions at small scales fl9| . 

To illuminate the multi-scaling at small scales, consider 
a picccwise continuous fracture surface with a number of 
vertical jumps of size e^, e.g. due to overhangs, micro- 
scopical defects or the grain size in granite and the fibre 
size in paper. On sufficiently small length scales S, the 
height variations in the surface will be negligible relative 
to the jump size t{. The structure function will therefore 
collect from each jump a contribution roughly propor- 
tional to bdl |2l|. Overall, the contribution to Eq. 
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FIG. 3: Statistical distributions of the height fluctuations 
P(A'h) sampled from a grid of 2000 lines in the direction per- 
pendicular to the fracture propagation in a 3D granite block. 
We show the distribution for logarithmically increasing length 
scales 8. Note that in addition, we have shifted the various 
distributions logarithmically for visual clarity. We plot on a 
semilog scale P(A' h)V2na^ versus (A'h) /\/2a 2 and observe 
at large scales a typical parabolic shape of a Gaussian distri- 

_ 2 

bution. The solid lines represent the curve y = e x and fit 
perfectly the experimental distributions above a characteristic 
length scale 8/S > S* ~ 50. Inset: The scaling behavior of 
the standard deviation of the distributions P(A'h(S)) allows 

us to extract the roughness of the fracture surface: a oc S 1 ' 



with C 



0.75. 



will be C k (S) ~ S 1 / k (Y Jl e k l ) 1/k where the sum is taken 
over all the jumps. We see that the structure function 
now scales with a /c-dependent Hurst exponent = 1/k. 
This behavior is observed for values of k close to or larger 
than unity. For small values of k the effect of the verti- 
cal jumps will diminish relative to the ordinary surface 
roughening and therefore £k ~ C f° r k *C 1. 

A few comments should be given on this /c-dependence. 
First of all, it is questionable to base arguments in favor 
of multi-affinity on the smallest scales. It is important to 
have in mind that the experimental scan of the surface 
at small scales might be spurious due to limitations of 
the profilomcters or scanners and contain artificial jumps 
due to overhangs etc. The discretization (coarse grain- 
ing) of the data also plays a crucial role, see [22]]. Fig. 
^demonstrates how sensitive the structure function is to 
the discretization. We generated graphs with roughness 
exponent £ = 0.7 from a fractional Brownian motion. We 
then filtered the graphs by representing the values of h us- 
ing 3, 5, 7, 8 and 9 bits (see Fig.^J. When decreasing the 
resolution (the number of bits), we observe at small scales 
a clear deviation from the expected scaling behavior and 
more importantly a separation of the structure functions 
Ck(8). The grain size in granite and the fiber size in the 
paper experiment may introduce a discreteness similar 
to that of sampling the fracture surface at a low resolu- 
tion. For example, in the granite experiment, we observe 
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FIG. 4: Influence of the coarse graining at small scales on the 
structure function C'k(8), (k=2,3,4), for synthetic fractional 
Brownian motions with a roughness exponent £ = 0.7. The 
100 samples have been coarse grained using different resolu- 
tions corresponding to represnting h by 3, 5, 7, 8 or 9 bits. 
This corresponds to using from r — 11 to 433 distinct values 
for spanning the range of /i-values. The grouping of the data 
sets into three different classes correspond to the three val- 
ues of k, with k = 4 at the top. Note that the separation of 
the structure functions on the small scales diminishes as the 
resolution is increased. 



a crossover from a microscopic multi-scaling to a pure 
sclf-afihic surface around the typical grain size (order of 
1 mm) , see Fig. [3J The multi-scaling, or more precisely, 
the separation of the structure functions at small scales, 
thus, most likely reflects the experimental limitations or 
the microscopic details of the material. 

In conclusion, we have analyzed experimental data on 
fracture profiles in widely different materials and have 
shown that the structure function ratios Rk (Eq. 0} con- 
verge to the values of a Gaussian process. We have ver- 
ified our findings by also computing directly the distri- 
bution of the height fluctuations and have shown that 
there exists a scaling exponent, £, permitting the rescal- 
ing of the height fluctuation distribution P(Ah(S)) ~ 
6~ <: G(Ah(S)/6 <: ). We find that the rescaling funct ion G 
has a Gaussian form. From a fundamental point of view, 
the distribution of the height fluctuations provides new 
important information about the morphology of fracture 
surfaces; information which is not covered by the calcu- 
lations of a roughness exponent. In addition, the wide- 
ranging convergence of the ratios Rk shows that there is 
no global multi-scaling of the structure functions. 
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